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Abstract 

We investigate the quantum motion of a neutral Dirac particle bouncing on a mirror in curved 
spacetime. We consider different geometries: Rindler, Kasner-Taub and Schwarzschild, and show 
how to solve the Dirac equation by using geometrical methods. We discuss, in a first-quantized 
framework, the implementation of appropriate boundary conditions. This leads us to consider 
a Robin boundary condition that gives the quantization of the energy, the existence of bound 
states and of critical heights at which the Dirac particle bounces, extending the well-known results 
established from the Schrodinger equation. We also allow for a nonminimal coupling to a weak 
magnetic field. The problem is solved in an analytical way on the Rindler spacetime. In the 
other cases, we compute the energy spectrum up to the first relativistic corrections, exhibiting 
the contributions brought by both the geometry and the spin. These calculations are done in 
two different ways. On the one hand, using a relativistic expansion and, on the other hand, with 
Foldy-Wouthuysen transformations. Contrary to what is sometimes claimed in the literature, both 
methods are in agreement, as expected. Finally, we make contact with the GRANIT experiment. 
Relativistic effects and effects that go beyond the equivalence principle escape the sensitivity of 
such an experiment. However, we show that the influence of a weak magnetic field could lead to 
observable phenomena. 

PACS numbers: 03.65. Ge, 04.80.Cc 
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I. INTRODUCTION 



Quantum mechanical systems in interaction with an external gravitational field are of 
particular interest, not only in theoretical physics, but also from an experimental point of 
view. Thirty years ago, the COW experiment |l| displayed the quantum mechanical phase 
shift experienced by neutrons waves due to their interaction with the Earth's gravitational 
field. The acceleration of a sodium atom due to gravity has also been measured, see Ref. 



^ound states of ultracold 
3j. Most of the theoretical 



More recently, the GRANIT experiment proved the existence o 
neutrons bouncing on a mirror in the Earth's gravitational field 
studies related to these experiments feature the Schrodinger equation in the linear Newtonian 
potential . This simple model gives satisfactory results in reproducing the experimental 
data because the particles or atoms used are nonrelativistic. 

The next step in the improvement of such models should be the inclusion of the spin 
degrees of freedom. In particular, the dynamics of neutrons should be described by the 
Dirac equation on curved spacetime. Many theoretical investigations have been devoted to 
the study of the spin-gravity couplings of a Dirac particle on the Rindler spacetime (see 
for example Refs. U El). F ew °f them focused on the problem of bound states of Dirac 
particles in a gravitational field (see, however, [8j]). The latter problem, which we investigate 



in the present work, is motivated by the GRANIT experiment but also has an intrinsic 
theoretical interest related to the understanding of gravitational effects in first-quantized 
systems beyond the equivalence principle j^. 

We focus on a neutral Dirac particle bouncing on a fixed mirror in a gravitational field. 
Three particular geometries are considered. The first one involves the Rindler metric, which 

n 

corresponds to the metric of a uniformly accelerated observer 10]. It is generally believed to 
describe a homogeneous gravitational field. The second metric with a cylindrical symm etry 
that we use is the Kasner-Taub metric, which is a particular case of Kasner metrics 
In contrast with the Rindler metric, it describes a genuinely curved spacetime, correspond- 
ing to the spacetime in the neighborhood of a plane brane jl^j], see also jl^j]. Because 
of the equivalence principle, both Rindler and Kasner-Taub metrics are equivalent in the 
Newtonian limit, but differ in relativistic corrections as we show. Finally, we also use the 
well-known Schwarzschild geometry. 

Our goal is to compute the energy spectrum of the Dirac particle. The energy is indeed 
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quantized because of the mirror which imposes particular boundary conditions. Using the 
Rindler metric, we can compute the energy spectrum in an analytical way. For the other 
geometries, we restrict ourselves to the first relativistic corrections that we compute by 
resorting to the usual relativistic expansion and by using a Foldy-Wouthuysen (FW) trans- 
formation It is of interest to compare these two different methods since the validity of 
the FW transformation in the present context is still subject to discussions, see Ref. pi and 
references therein. We also allow for the neutral Dirac particle to be nonminimally coupled 
to a weak magnetic field and show that some observable effects could be derived in future 
experiments. 

Our paper is organized as follows. In Section |H] we introduce the notation, write the 
Dirac equation on curved spacetime and focus on cylindrically-symmetric geometries. We 
also discuss the boundary conditions encoding at best the presence of the mirror on which 
the particle bounces. In Section IIHI we proceed to the same analysis with a spherically- 
symmetric geometry. In Section IIVI we particularize the metrics that we use: Rindler and 
Kasner-Taub. Then, we solve the Dirac equation on those spacetimes. On the Rindler 
spacetime, the Dirac equation can be solved analytically, taking into account the boundary 
condition induced by the mirror. Moreover, we compute the energy spectrum of the particle 
in the Rindler and Kasner-Taub spacetimes using a relativistic expansion. In Section 
the same computations are done, but this time by resorting to a FW transformation and 
also considering the Schwarzschild geometry. Finally, we apply our results to neutrons in 
the Earth's gravitational field and compare them to the GRANIT experiment in Section fVTl 
Our conclusions are given in Section I VI II 



II. DIRAC EQUATION ON A CYLINDRICALLY-SYMMETRIC SPACETIME 
A. Dirac equation in curved spacetime 

We introduce the notation and recall the main theoretical features that we will need 
about the Dirac equation on curved spacetime. First of all, the Clifford algebra is given by 
{7a, 7b} = ^Vab where a, b = 0, . . . , 3 . The metric is given by ds 2 = r] a b 9 a <g> 9 b , where 9 a are 
the coframe fields that diagonalize the metric. We adopt the "mostly minus" convention 
r] = diag(H ) and use units such that h = 1 = c. 
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The covariant derivative for a two-spinor in a curved background has been known for a 
long time ^(|. When acting on a Dirac spinor it is given by 

V¥ = dtt - i U ab Z ab V = (d + , Zab=l [la, lb] , 

where H ab are the generators of Lorentz transformations and uj ab are the usual Levi-Civita 
connection one-forms. An important property of V is that it preserves the Dirac matrices, 
i.e. 

V7 6 = dj b - u c b 7 C + [Q, 7 6 ] = . 

Finally, the Dirac equation on a curved background reads 

i7 a e A * a V M * ~ rn^ = z7 a V a * - = , (1) 

where the symbols e M a denote the coefficients of the vierbeins e a in coordinates. In the 
following, Latin indices will denote orthonormal indices, whereas Greek indices will denote 
curved (holonomic) indices. Hated specific indices will also be used to denote flat (non- 
holonomic) indices. 



B. Explicit form 

The metric 



ds 2 



A 2 {()dt 2 - B 2 {Q{dx 2 + dy 2 ) - dC 



(2) 



clearly describes a spacetime with cylindrical symmetry along the ( axis. For the moment, 
the functions A and B are arbitrary functions of ( ; they will be specified later on. Knowing 
the metric, it is straightforward to write the (co)frame fields: 







9° = 


Adt 




< e 3 = , 


< 


03 = 


B dx^ , 


J = 1,2 


k e 3 = d ( 




6* = 


d( 





From the zero-torsion condition, we deduce the only nonzero components of the spin con- 
nection one- form uA : 



the prime denoting a derivative with respect to the variable ( . 

We want to study a neutral fermion of spin 1/2 with mass m and anomalous magnetic 
moment /i n , not only in a gravitational field, but also in a weak external electromagnetic 
field described by the field strength F^ v . Including such an electromagnetic interaction is 
achieved by adding to the Dirac equation a nonminimal coupling term of the form —fi n Il a bF ab 



17| . If the electromagnetic field is purely magnetic, this term reduces to {i n T,^B k , with B k 
the magnetic field {k = 1, 2, 3) and S fc the spin matrices. Using the previous results, we 
have explicitly from Eq. ([T|) 

r - 7 - A' 2B' -i 

i-fA-% + i^B- 1 ^- + ij 3 d c + - 7 3 (-r + — ) - m + fi n Z~ k B k * = . (3) 

Z A tS 

Concerning the metric (J2J), it is worth mentioning that we will consider it independent of the 
magnetic field. The problem of writing cylindrically symmetric solutions of the Einstein- 
Maxwell equations is a nontrivial one, out of the scope of the present work. We assume here 
that both the gravitational and magnetic fields are weak enough so that one can neglect 
their mutual interactions. 

C. Dirac wave function 

We are interested in solutions described by the following positive-energy ansatz: 

^ = e -itot e ik-x A ~i/2 B -i e ^ u > . 

Equation (jSJ) becomes 

u^A- l e - ftB- 1 ® + rfe' - me + ^ k B k e = o (4) 

where k = (k 1 , k 2 , 0) and $ = k^j with j = 1,2. 
We now define the spinor U a such that 

UlU a = l, 7°U a = U a , h*U a = i<TU a , (5) 

with f, = ft/k, k = \k\ , and a = ±1 . In general, U a defined by the relations (0) is not an 
eigenvector of the magnetic term S^.i3 fc . It will be such if 

~h\ = , 
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or equivalently if B is orthogonal to k, with B 3 = . Thus, if k = (k, 0, 0), we will assume 
that B = (0, B, 0). Then, the following relations also hold 

ftU a = ikaj 3 U a , 
2 

As U^ 3 U = , it is convenient to decompose 6 as follows: 

B = F(C)U + G{crtu, (6) 

where we dropped the dependence in a in order to simplify the notation. The normalization 
condition 

y"e+edc = i (7) 

implies that the following relation must hold: 

J (\F\ 2 + \G\ 2 ) d( = 1 . (8) 

Moreover, the application of the condition (jZJ to Eq. (pij allows us to write the energy of 
the neutron as 

Q^H^QdC, (9) 



with 

H = —AB~ l f, + iA^d c + mA t + // n S^#. 

In principle, Eqs. (jHJ), © and a boundary condition (discussed below) are sufficient to 
compute the energy spectrum, provided we know the functions F and G . We thus have to 
find the differential equations that these functions satisfy. 

Due to the decomposition @ and to the fact that the spinors U and 7 3 [/ are linearly 
independent, Eq. (J3J splits into two parts: 

D-F — i<pG — iG = 0, 

(10) 

-D + G -icpF + iF' = 0, 



with 



D± = ^- 1 ±mT^ I V = B- 1 ka. (11) 
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Once F is known, G is readily given by 



G = —{d c -<p)F. 
From Eqs. (jlOjl. one can extract the following equation 



F "-it F ' + 



D + D_-<p 2 + ^<p-<p' 



F = 0. 



Finally, upon rescaling F according to 



F = Df /(C) 



we are left with 



D' 



D + D_ -<p 2 + -±(p -(f/ + 



D" 3D 



12 



f = o. 



(12) 



(13) 



(14) 



(15) 



D + T r ' 2D + W\ 

Together with Eqs. (fT2"|) and (jTljl . the above equation enables us to find F and G once / is 
known. 



D. Boundary condition 

The particle we are considering bounces on a mirror placed at £ = . This situa- 
tion should provide us with a particular boundary condition and should eventually lead to 
the quantization of the energy. By analogy with the standard Schrodinger equation, one 
could demand that the wave function vanish at the mirror. This condition can easily be 
implemented with a Klein-Gordon field 181 and indeed leads to the existence of bound 

n 

states |19|]. However, for a Dirac particle, things are not so simple. For example, imposing 



01 



C=o 



is not very fruitful since the relation fl2|) would lead to the conclusion that G 



is zero everywhere. The same conclusion holds if one requires that the probability density 
jo _ QtQ _ |_p|2 _|_ |£»|2 vanish at the boundary. Moreover, the j 3 -component of the current 



is given by 



f = 7 3 = —(FF 1 * - F*F') 



Since F and F' are determined only up to an arbitrary complex coefficient and thus can be 
chosen to be real functions (see Eq. j 3 is identically zero. Physically, this expresses 

the fact that a repeatedly bouncing particle is described by a stationary state. 
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A satisfactory first-quantized boundary condition can be found if we represent the mirror 
by a scalar potential 20j] V such that V|^ = V ^> m and ^|^ >0 = . Formally, at the level 
of the Lagrangian, the introduction of such a potential amounts to replacing m by m + V . 
Then, in the region ( < where the potential blows up, Eq. (|15j) becomes approximately 

f" « V 2 f , 

which yields 

F^N Vo /2 e VoC , G « iN V? /2 e VoC , (16) 

where iV is a normalization constant. Since, at £ = 0, the wave function has to be continuous, 
we assume from Eq. (116)1 the Robin boundary condition 

^l c =o = ~i G\c=o (17) 
which, using Eq. (fT2|l . can be rewritten as 

(D + + ^)F| f=0 =F'| c=0 . (18) 

Thanks to this last relation, the energy of the neutron will be quantized. The presence of 
ip ensures that only two solutions exist: The energy will be spin-dependent. Moreover, at 



leading order, Eq. (fTSjl becomes F\^ =Q w . This is the boundary condition we are used to 
see when the Schrodinger and/or the Klein-Gordon equations are employed |igL I2I | . 



III. DIRAC EQUATION ON A SPHERICALLY-SYMMETRIC SPACETIME 

We now turn our attention to spherically-symmetric metrics. Such a metric reads 

ds 2 = e 2 ^dt 2 - e 2A ^dr 2 - r 2 du 2 (19) 
and the corresponding Dirac equation is 



\ 2 r J r \ 2 J rsmd 



^ = 0, (20) 



where a prime now denotes a derivative with respect to the radial coordinate r . It is 
convenient to use the following representation of the Dirac matrices 



7° 






Decomposing the Dirac spinor into a pair of bi-spinors 

* = 

we find that Eq. (|20j) yields the following two equations 

i e-*d tX + (T 3 e~ A d r r] + -$r) + ^ A ( $' + - J a 3 rj = m X , (21) 

r 2 \ r J 

-i e~^d t rj + a 3 e~ A d rX + -$x + ^ e " A ( & + -) a 3 X = mr] , (22) 

r 2 \ r J 

where Tp is the Dirac operator on the unit sphere, 

$ = a 1 de + -^—;<J 2 d v + \cot6a 1 . (23) 
sin 6 2 

In order to solve these equations, we may use the well-known decomposition of bi-spinors 
in terms of spherical spinorial harmonics Q, presented from a group-theoretical point of view 



in (almost) all text books on relativistic quantum mechanics (see for example Ref. 22 1). 
Nevertheless, for pedagogical purposes, we find it useful to reconsider this decomposition 
from a more geometrical point of view, the method having the advantage of being easily 
extended to more general situations (higher dimensions, other geometries, etc.). Let us 
separate the variables and write the bi-spinors x an d f as products of the positive-energy 
exponential e~ tuJt multiplied by radial functions and spherical spinorial harmonics: 

x = e -^F(rM6, if), V = -ie- <wt G?(r)fi(0, <p), (24) 

where Q has to be specified appropriately. By inserting this ansatz (|24|) into Eqs. ()21|) and 
(J22J), it is easy to see that the Dirac equation is satisfied provided 

(i) we define Q = a 3 Q ; 

(ii) the spinor Q verifies the equation 

a 3 lpn = Ktt. (25) 

By integrating fl^a 3 ififl by parts on the sphere, it can be found that the eigenvalues n are 
real numbers. On the other hand, we know that the eigenvalues /i of the Dirac operator ()23j) . 
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are purely imaginary. From an eigenspinor of the Dirac operator we obtain a solution of 
Eq. (|25|) by taking the combination Q = T + ia 3 T, the corresponding eigenvalues being 
related by k = — i fi . 

For completeness, let us recall that the eigenspinors of the Dirac operator (on any sphere 
of arbitrary dimension) can be obtained by using a Killing spinor and the usual spherical 
scalar harmonics Y as follows. A Killing spinor K is a spinor verifying the equations 

D & K k = ±-a & K k . 

On S 2 there exist four Killing spinors, k = 1,...,4 j^J. Taking anyone of them, 
we may obtain the eigenspinors of the Dirac operator by considering the combination 

= daX a a K k + XY K k [24]. Substituting this ansatz into the Dirac equation and using 
the well-known spectrum of the spherical scalar harmonics, we obtain A = i fi with fi G Z . 

The well-known equations that the radial functions F(r) and G(r) have to satisfy are 

£>_F -%G' -iyJG = 0, 

(26) 

-V + G + iF' + i V+ F = 0, 
where 

^ A— $ _i_ A $' 1 ± ne A 

V± = uj e ± me , (p± = — + 



2 r 

These equations can be seen as the counterpart of Eqs. (fHH) . From Eqs. (J2HI), we obtain that 
F and G are given by 



(27) 



G = ^-(d r + v + )F, (28) 

these equations being similar to Eqs. (j!3p and (|T2*|) obtained previously in the cylindrically- 
symmetric case. 

The number n is a nonzero integer which can directly be interpreted in the Schwarzschild 
geometry. It represents the eigenvalues of the operator constructed on the time translation 
Killing vector and on the Penrose-Floyd tensor j^j, see also Eq. (23) of Ref. j^ . 

Due to the symmetry of the problem, the boundary condition (j!7|) will be slightly modi- 
fied. Let us indeed consider that the scalar potential is equal to a large constant Vq in the 
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region r < R, with R an arbitrary radius, and vanishes for r > R. Then, we find from 
Eqs. (JHJ) that 

F\ r=R =-ie A G\ r=R . (29) 
In the Schwarzschild geometry for example, we have 

F\ r =R = -^T=^ G l=R ■ (3°) 
V ~R 

Note that, obviously, such a boundary condition is only valid in the region R 3> r s . We 
stress that the present first-quantized formalism is no longer valid near the Schwarzschild 
radius, where a second-quantized formalism is definitely needed. 



IV. EXACT SOLUTIONS 

A. Geometrical preliminaries 

Before proceeding to the resolution of our equations, we have to specify the metrics we 
will use. Usually, the flat Rindler metric is used to describe the physics of a uniformly 
accelerated system in flat spacetime (see e.g. [27] chapter 6) 

ds 2 = (1 + g() 2 dt 2 - dx 2 - dy 2 - d( 2 . (31) 

With the notation introduced in Eq. (J2J), it corresponds to 

A=(l + gQ, B = l. (32) 

Another spacetime metric, that describes a cylindrically symmetric solution of the vacuum 
Einstein equations, is given by the Kasner-Taub metric [l2| 

ds 2 = (z/z )- 2/3 dt 2 - (z/z ) 4/3 (dx 2 + dy 2 ) - dz 2 . (33) 

It can be interpreted as the gravitational field around an infinite plane . The constant 
z can be related to the gravitational acceleration in the Newtonian limit. In this limit, 
we set z = Zq — ( and expand the metric component goo to first order in (/zq, which gives 
goo = 1 + |f- = 1 + 20jv • The Newtonian potential is thus 4>n = g( where g = | z . 
Considering the full relativistic metric ()33|) 

ds 2 = (l-3g Cr 2/3 dt 2 - (1 - 3 g C) 4/3 (dx 2 + dy 2 ) - dC, 2 (34) 
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corresponds to taking 

A=(l-3sCr 1/3 , £ = (1 - 3sC) 2/3 (35) 

in Eq. (J2J). In order to restore the c factors, we have to replace g by g/c 2 and t by ct as 
usual. Let us note that the metrics defined through Eqs. (j3*2~|l and (jHEJ) differ at order c~ 2 . 
For the sake of simplicity, we shall discuss both simultaneously using 

A= 1 + 4 , B = l-b^, with a =-1,1, b = 2,0, (36) 

for the Kasner-Taub and Rindler metrics, respectively. 

B. Bound states in Rindler metric 

Exact solutions of the free Dirac equation on Rindler spacetime are well-known (see for 
example Ref. [3], or sometimes in a disguised form, Ref. [2^ and Refs. therein). In the 
present section we particularize these solutions to our problem, namely a particle bouncing 
on a mirror, and show that boundary condition (j!7)l implies the quantization of energy. 

Instead of the spinor U a of Section III CI we introduce the spinor W a obeying 

j l ^ 2 W a = iaW a . 
The positive-energy Dirac spinor is now defined as 

^ = e — V^e CT (p) 

where p = ( + - and 

Q a (p) = \-(p)W a + \ + (p) 1 d W a . 
It can be shown that the functions A ± obey the following equations [28] 

[pd p pd p ] A ± -- 



1 \ 2 
too . 1 



(m 2 + A;V+(-± 9 
whose solutions are 

A± = H il ) / g ±l/2^ K P) > 
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A*, 



where is the Hankel function of the first kind (the Hankel function that vanishes when 
p is going to infinity) and k 2 = m 2 + k 2 . 

The spinor is an eigenspinor of the operator 7 X 7 2 . Since this operator commutes with 
the Hamiltonian, o is a good quantum number. Depending on its values, the large and small 
components of CT can be identified with the functions F and G by the relations 

FocA~ + aA + , GocA--aA+. 

The position of the mirror being given by p — 1/g and because A - = (A + )*, the boundary 
condition (j!7|) can be rewritten 

m£ /9+1/2 (iK/g)} + o- m2 /9+1/2 (™/g)} = o . (37) 

Let us emphasize that if the spin contributions (the 1/2 in the index of the Hankel functions) 
are neglected, these Hankel functions become real and Eq. (|37j) yields 



H 



(i) 



> /g (iK/g) = o 



(38) 



which is the usual boundary condition in the Klein-Gordon case A plot of the left-hand 
side of the boundary conditions (|37|) compared to (|3*Hj) is drawn in Fig. Q where we defined 
e = uj — k c, restoring the c factors. Obviously, it leads to positive, quantized energy levels 
£ nu . In the nonrelativistic limit, these levels all coincide. This is illustrated in Fig. El 

In order to discuss this nonrelativistic limit, it is useful to introduce the following quan- 
tities 

UJC 1 IK C 2 

u — i h - , w — . 

g 2 gu 

For weak gravitational fields, we have \u\ ^> 1 and consequently \w\ « 1 . In this case, we 
can assume that e/kcC 1 and approximate the Hankel function by [30! - Eq. (9.3.37)] 



H$p(uw) oc Ai 



mg z 



1/3 



(_9_Y /3 



+ ~ { T- ) Ai 

c \4m/ 



mg' 



1/3 



+ 0(c 



(39) 



where Ai denotes the regular Airy function. Note that KRimcin the nonrelativtistic limit. 
Moreover, the imaginary part of can be neglected in this limit, so that the condition ()37|) 
together with Eq. ((HH)) gives 

„2\ 1/3 

a n , (40) 



mg' 
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with a n the n th zero of the regular Airy function. These can be found for example in Ref. [30 



Table 10.13]. A WKB approximation of these zeros can also be found in [3y, Eq. (10.4.94)]: 



a, 



y(n-l/4) 



2/3 

(41) 



This approximation is precise up to 8%. Note that the asymptotic values of the energy levels 
e na in Fig. El are given by —a>i ~ 2.34 and — a 2 ~ 4.09 as expected from Eq. P0|). 

C. Relativistic expansion 

1. Lowest order solutions 

Although it is possible to find an analytical solution of the Dirac equation in Rindler 
spacetime, the determination of the energy is rather problematic in the other cases, due to 
the nature of the singularity in the differential equations. To our knowledge, no explicit 
solution in Kasner-Taub geometry has been found yet. It is thus interesting to expand the 
equations in powers of 1/c 2 . Moreover, this procedure will allow us to find a solution in 
Kasner-Taub spacetime and to study the influence of a weak magnetic field. The equation 
we need to solve, that is Eq. ([To]), can be expanded up to the order c~ 2 with the appropriate 
choice of A and B given by Eqs. (JHHJ). We define the energy uj as 

k 2 

uj = mc 2 + +g. (42) 

2m 

^From Eq. (|15j) . we obtain at lowest order 

-& + m9Cf ~^ Baf = £L (43) 

which is the expected Pauli equation with a linear gravitational potential and a magnetic 
field. To go ahead, we simply assume that the magnetic field is constant in the y direction, 
i.e. 

B = B . (44) 

This choice is relevant in the context of the GRANIT experiment [3| since it can easily 
be added to the current experimental setup. Equation (|4*5|l is then simply a Schrodinger 
equation with a linear potential, with the boundary condition (fT%|) given at this order by 
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f\(=o = • The solution of sucn an eigenequation is well-known |2l|, Problem 40] 



£na = - (^j- ) a n -^-B a, (45a) 

= ,, Ai [9 C + «„] , (45b) 

|Ai(a n )| 

with 

6=(2m 2 g) 1/3 . (45c) 

Without magnetic field, Formula ()45a)l tells us that the energy of the particle, and thus 
the height at which it can bounce, only depends on its mass and of the strength of the 
gravitational field. If a nonzero magnetic field is present, the energy depends on a. As it 
is deduced from Eqs. ©, cr/2 can be interpreted as the spin of the particle along the y 
direction. This will be detailed in the following. 

As a consistency check, we note that the formula (jHJ), taken in the nonrelativistic limit 
where F = f and G = , gives 

The average values are computed with respect to the wave function (|45b|) . The expressions 
for these average values are analytical, as it can be seen in Eqs. PH]) . 

2. Relativistic corrections 

We are mainly interested here in the relativistic corrections to the energy spectrum. These 
can be obtained by expanding the formula (jHJ) at the order c~ 2 . The wave functions which 
have to be used are the properly normalized F and G computed thanks to Relations ©, 
(H2J) and (JT3J), with / given by (|45b|) . We can keep the same boundary conditions as in 
previous section. After some algebra, we find 

uj = ujq + uj 2 + 0(c -4 ) , 

with ojq the nonrelativistic energy pE|l and 



(k 2 + &))'' 



U2 = i - ^ — + 



M> + (1+2()) (o^ 

2m v J 2m 



I + (l-.)!^ + (l + Ji )5. (47) 
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The various relativistic corrections appearing in the above equation are now interpreted in 
analogy with the energy E of a classical particle on a curved spacetime. With the metric (J2J), 
we obtain 



and thus 



E = AJm 2 c 4 + B~ 2 c 2 k 2 + c 2 p^ 



2 



Pc 

E m mc H 1 — + mgQ 

2m 2m 

8m¥ + c'2m +l + V 2m + 1 ' c 2 ' 1 J 

Apart from the spin (a) dependent term, Eq. fljTj) clearly reproduces the terms appearing 
in the classical formula (|48p. The first term in co> 2 is the usual relativistic correction arising 
from the expansion of the relativistic kinetic energy in powers of 1/c 2 . The next two terms 
between square brackets are redshift corrections. The first one gives a quantum mechanical 
translation of the usual redshift formula for a particle falling in a gravitational field, viz 

E{L) = (l + ^) E(0) . 

The second is an additional transverse redshift term. The ( 2 term is a curvature correction. 
As expected, it vanishes in the flat Rindler geometry. The last term is a spin- dependent 
correction. 

The properties of the regular Airy function lead to the relations 

<?<} = -o 3 (0 - »V„ (C?c> = -o 3 (O - 2q » (0 . ( 49b ) 

which allow to express Eq. (}4Tj) in terms of a n (}4*Tj) and 6 (145 c|) . 

V. FOLDY-WOUTHUYSEN TRANSFORMATION 
A. Kasner Taub and Rindler metrics 

It is interesting to compare the previous energy spectra with the results obtained by 



resorting to a FW transformation }l5 . 
dealing with is indeed still subject to discussions 

17 



31| . Whether this procedure works in the case we are 



The Dirac equation Q can be recast in the Hamiltonian form 



(50) 



where 



H = (5m- fi n (3Z k B k + + P, 
O = AB-^pj + Ac?p c - % - (a' 



(51) 



2B'A 
B 



a 



1,2, 



V = m(A-l)p. 
We recall that the operators O and V satisfy 

{P,O} = 0=\/3,V] 

and that pk = —idk ■ Using the standard Bjorken and JDrell conventions, at the first order 
in 1/m, the FW Hamiltonian computed from (JoTj) is 

1 



O 2 

H FW = (3m + (3— + V 
2m 



[O, [0,V]} 



The positive-energy part of this FW Hamiltonian is given by 

k 2 



H 



FW 



mc 2 + 



1 + (1 + 26) 



g ■ x 



2m 



+ 



1 + 



g ■ x 



2m 



+ mg ■ x — jjL n S ■ B 



g-p 



(52) 



2mc 2 c 2 2mc 2 ' 

with g = (0,0, g). The lowest order FW Hamiltonian is almost identical to the one which 
can be read from Eq. (|4H|) provided that one sets S y = cr/2 . This confirms the interpretation 
of a/2 given in the previous sections. Moreover, the relativistic corrections are on average 
identical to those obtained in the formula (|4"7)l since (p^) = 0. As the FW transformation is 
only performed up to the order m _1 , we miss the kinetic corrections in m~ 3 . So, comparing 
the method used in the previous section with the FW method at the same order in 1/m, we 
may conclude that both approaches lead to the same results. Let us also note that for the 

n 

Rindler metric without magnetic field, we recover the result of Ref. p. 



B. Schwarzschild metric 



Since we have shown that the usual relativistic expansion agrees with the FW technique 
for Kasner-Taub and Rindler geometries, we apply the latter to the Schwarzschild geometry. 
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The Hamiltonian corresponding to Eq. (|20|) is 



H = -ia e 



3„$-A 



$' 1\ *e* ( „ cotO 

I — zo;^ 

r r I r \ 2/ r sin 



d r + — + -)- ia 1 — d e + 



ds + me®/3 . 



After a FW transformation, we find that the positive-energy part of the Hamiltonian is 



if F \v = m c 2 



3g-r 



Pr 

2m 



9 ■ r 



2mr 2 



+ mg ■ r 



S 



+ i 



(53) 



2mc 2 ' ' ' 2 mc* 

where g = —GMl r /r 2 and £(£ + 1) is the eigenvalue of the squared orbital momentum. 
Formula ([53)1 agrees with the corresponding formula (28) of Ref. j^l upon performing a 
change of coordinates from isotropic to the usual Schwarzschild coordinates. Here again, 
(p r ) = and the last term in Hamiltonian ()53|) brings no contribution. 

Let us point out that the classical energy of a nonrelativistic particle moving in the 
spacetime (|T9*jl geometry is given by 



E — e*\/ m 2 c A + e 2A c 2 p 2 + c 2 — 



In the particular case of the Schwarzschild metric, we have 



E pa mc + 



Pr 



+ 



2m 2mr 2 



+ mg ■ r 



(p 2 + L 2 /r 2 ) 2 „g-rp 2 g-f L 2 m (g ■ f) 2 



(54) 



8m 3 c 2 c 2 2m c 2 2mr 2 2 c 2 

The classical energy given by Eq. (|54jl is again identical to the Hamiltonian (jH^j) apart from 
the spin-dependent terms and the kinetic corrections in m~ 3 . 



VI. COMPARISON WITH THE GRANIT EXPERIMENT 

What can be measured experimentally is the critical height, corresponding to the classical 
turning point of the neutron. It is given by 



£na = mgh r , 



(55) 



which can be rewritten in our approximations 



6 2mg mg 
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(56) 



At the lowest order, and without magnetic field, we simply obtain 



(57) 



This relation has been successfully checked for the first two bound states by the GRANIT ex- 
periment p. Indeed, knowing that 6 = 0.17/im -1 for a neutron in the Earth's gravitational 
field, the predicted heights are 

hi = 13.7 /im, /i 2 = 24.0/im, (58) 

while the experimental results concerning these states are 

h? p = 12.2 /im ± 1.8 syst ± 0.7 stat , hf = 21.6 /im ± 2.2 syst ± 0.7 stat . (59) 

If we still stay at the lowest order but allow for a weak magnetic field, then, formula (JHHJ) 
predicts a splitting of the critical heights depending on the values of a. A given level h n shall 
be split in /z n l and h n _i , both levels being separated by a value which does not depend on 
n, that is 



A/i« . (60) 

I m 9 

In the future of the GRANIT experiment, it seems possible to deal with a magnetic field 
whose density v is of order O,r/ m 0. As th e eh— c size of the experiment is 
around 10 /im, a value such as Bq ~ 10~ 6 T could be produced (we assume the experimental 
apparatus to be shielded from the Earth's magnetic field). In such circumstances, one would 
observe Ah ~ 0.6 /im. An experimental confirmation of this point requires an increase of 
the experimental accuracy, which should be possible in a near future [32] . Note that a value 
such as v = 0.1 T/m is still in the weak magnetic field regime since ii n v/mg ~ 6 %. 

The relativistic corrections are particularly interesting, since they involve couplings be- 
tween spin, momentum, and gravity. In particular, there exists the term S ■ (g x p), which 
states that even without magnetic field, the particle will bounce at different heights depend- 
ing on the value of its spin. Again, this term will split a given height-level into two separate 
levels, distant by the quantity 

Ah sg oc ^ , (61) 

where Vh is the horizontal speed of the particle. In the GRANIT experiment, v\ t ~ 6.5 m/s jsj. 
This leads to Ah sg = O(10 _18 /im): The experimental detection of this phenomenon unfor- 
tunately seems hopeless with such an experiment, as well as the detection of the other 
relativistic corrections, as pointed out in Ref. 
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VII. CONCLUSION 



In this work, we obtained the quantization of the energy of a Dirac particle bouncing 
on a mirror on a curved spacetime. Consequently, the critical heights at which the particle 
bounces are also quantized. This is due to the presence of the mirror, which imposes a Robin 
boundary condition leading to the existence of a discrete set of bound states. We computed 
the spectrum of these bound states for several choices of spacetime metrics. Firstly, we solved 
the problem analytically on the Rindler spacetime. We then computed the energy spectrum 
in Kasner-Taub and Rindler geometries at the second order in 1/c 2 thanks to a relativistic 
expansion. At lowest order, these energies are equal in both geometries - in agreement with 
the equivalence principle. However, they differ in the relativistic corrections through redshift, 
curvature and spin terms. One of the most conceptually interesting corrections is a term in 
S • (gx p), which shows that the energies, and thus the critical heights, are spin-dependent. 
We also expressed the Hamiltonian in Rindler, Kasner-Taub, and Schwarzschild geometries 
using a Foldy-Wouthuysen transformation. We checked that the results are identical for the 
first two geometries, showing a posteriori that the Foldy-Wouthuysen transformation is valid 
in this framework. Finally, we compared our results with those of the GRANIT experiment. 
The nonrelativistic approximation is enough the reproduce the current experimental results, 
but the relativistic corrections appear too be small to be detected. Moreover, we showed 
that the presence of a weak constant magnetic field leads to observable effects. 
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FIG. 1: Numerical plot of the left-hand side of the Dirac boundary condition Q37|) for a = +1 (solid 
line) and a = — 1 (dotted line), and of the Klein-Gordon boundary condition (|38jl (dashed line). 
The zeros of these curves (full circles) are the energy eigenvalues, e na , given in units of (ng 2 /2c) 1//3 . 
We fixed k = 2, g = 1 and c = 1. 
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FIG. 2: Numerical plot of e na in units of (n g 2 /2c) 1 ^ 3 versus the dimensionless parameter (k c 2 /g) 1 ^ 
for the first two energy levels. We fixed g = 1. 
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